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ABSTRACT 


iIntensatiy Scimtillatien im a laser beam at 0.63 
micrometers in the marine boundary layer has been studied 
over a 4.3 kilometer horizontal path across Monterey Bay 
and also from shore to ship at San Nicolas Island. Optical 
and micrometeorological measurements of the refractive 
iImedex structure constant, Cy agree to approximately one 


standard deviation giving values in the range 2.0 x jones 


-_ ¥ = , 
moms. 0x 10° (eters) ase Optical and meteorological 
data.correlation improved as mean wind speed increased. 
Evidence was found to support the reported existence of a 


Saturation region for the laser beam logarithmic amplitude 


Variance. 
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Tice OR ODUCT TON 


Atmospheric effects on laser beam propagation are 
divided into three general categories: refraction, absorp- 
tion, and particulate scattering. This investigation was 
exclusively concerned with atmospheric refraction effects. 
Refraction may be further subdivided into large-scale 
refraction, where the beam's spatial path is altered, and 
small-scale refrvyaction, where the beam is differentially 
focused and de focused. Large-scale Tetrection, commonly 
referred to as beam wander, was found to be negligible. All 
research efforts were directed toward the small-scale refrac- 
tion phenomenon commonly referred to as beam or atmospheric 
scintillation. 

Bett) OF purpose of the research was to make and to 
correlate direct optical measurements of laser beam scintil- 
lation with meteorological measurements of certain atmospheric 
parameters. Using Rytov's method as described by Tatarski 
[1], the optical and meteorological measurements of beam 
scintillation were compared by calculating the atmospheric 
refractive index structure constant which statistically 
describes atmospheric refractive index fluctuations. The 
research team attempted to predict laser beam scintillation 
by making micrometeorological measurements of atmospheric 
parameters. 

me proyect effort consiSted of two primary areas, an 


optical determination of the refractive index structure 
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constant (Cr) and a meteorological determination of the 
refractive index structure constant (Cre) - These measure- 
ments were made simultaneously and compared for correlation. 
iaiseexperiment differs from previous experiments conducted 
by Livingston [2], Ochs [3] and Fitzmaurice [4] in that the 
laser path was over the ocean instead of over land. 

This research project was funded by the Naval Ordnance 
Laboratory and was divided into two related areas, meteorol- 
ogy and optics. |The meteorological research group, headed 
Dy Dr. Komiewe fetadsen of the Naval Postgraduate School, 
made and evaluated the micrometeorological measurements of 
memospheric turbulence. The optical research group, 

Meaded by Dr. Eugene Crittendon (until his departure for 
Australia in January) and later by Dr. Alfred Cooper, made 
and evaluated the optical measurements of laser beam 


Seancillation. 


= 
= ’ = on % 
Cae SO ong Soe 
= * My an 
ate ce 


pwe to the wide extent of exeiwor the areas of investiga- 
tion, each group worked independently of the other to obtain 
an estimate of. atmospheric turbulence. These estimates were 
compared for correlation after data reduction and analysis. 
The inherent nature of this organization promoted speciali- 
zation of the personnel involved in either the optical or 
meteorological areas of research. The author, by virtue of 
his background in electro-optics, concentrated ils Sabie ye cs 
in optical research. Hence, this report is devoted pri- 
marily to the procedures and conclusions of the optical 


research group. This thesis represents a comprehensive 
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Mescription of the efforts and results of the optical 
propagation study. Meteorological techniques and experi- 
mental results are described briefly when they relate to 


optical data. 
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II. PROJECT BACKGROUND 


A. HISTORICAL SUMMARY 

The earliest observations of scintillation in a light 
beam propagating through the atmosphere were made by ancient 
astronomers who noticed the "twinkling" of stellar images. 
Stellar scintillation was measured quantitatively by Middleton 
mrer949 [5] and 1952 [6]. Chandrasekhar, in 1952 [7], 
estimated the phase and amplitude fluctuations of starlight 
caused by atrospheric turbulence by means of a geometrical 
optics approach. In 1953, Mintzer [8] applied the Born 
approximation to the scalar wave equation to describe the 
propagation of sound waves through a turbulent ocean, a 
problem similar to the propagation of light through a turbu- 
lent atmosphere. 

Mennore and Wheelon [9], in 2955, used Chandrasekhar's 
approach to study line of sight propagation of radio waves 
mecne troposphere. They calculated the variance and 
Spatial correlation functions for the amplitude and the 
Bilis e of the propagating wave. Two articles by Wheelon 
published in 1959 [10, 11] review the state of the art of 
Wave propagation at that time. These articles cover 
turbulence theory and atmospheric propagation of radio 
waves which is a similar problem to propagation of optical 
beams. 

In 1960 and 1961, the Russian monographs by Chernov [12] 


and Tatarski [1] were translated into English. These works 
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described the method proposed by Rytov [14] (see Section 
mii.) for the derivation) for analyzing the diffraction of 
light by ultrasonic waves ‘and later applied by Obukov [15] 
to the propagation of light waves through a turbulent 
atmosphere. Strohbehn [16] presents a good discussion on 
the range of validity and approximations made by Rytov's 
method. The calculations in this project were based on 
Tatarski's formulation of Rytov's method. Rytov's method 
takes into consideration multiple scattering and diffrac- 
muon efiects Pha —e represents a significant advance- 
ment in the theory of atmospheric propagation of planar or 
@emerical electromagnetic waves. Strohbehn's article, 
meshed in 1968, is an excellent review of atmospheric 
propagation theory. 

Lumley and Panofsky [17] wrote (in 1964) a monograph 
covering many facets of turbulence theory as applied spe- 
cifically to the atmosphere. This defined the meteorological 
theory used by the Naval Postgraduate School's research group 
to obtain estimates of the atmospheric refractive index struc- 
ture constant. In 1972, Davidson, in a Naval Postgraduate 
School report [18], comprehensively described the meteoro- 
logical measurements made and the theory applied in this 
project. Articles by Wyngaard (et al.) in 1971 and Panofsky [19, 
20] in 1968 provide further references on meteorological 


turbulence theory as applied by the research group. 
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me RECENT DEVELOPMENTS 

In his monograph, Wave Propagation in a Turbulent 
Medium [1], Tatarski reported the results of some of the 
earliest experiments that compared micrometeorological 
estimates to optical estimates of atmospheric turbulence 
Using Rytov's method. These experiments were performed in 
the late 1950s using an incoherent light source. The 
propagation path ranged from 250 to 2000 meters over a 
Missian steppe. “Fxperimental results strongly supported 
his Pemmiat.on of propagation theory. The fluctuations of 
beam intensity and phase caused by atmospheric turbulence 
were found to have a log-normal probability distribution. 
Tatarski's description of the scintillation frequency spectrum 
and the dependence of the fluctuation variance on path length 
were found to match propagation measurements. 

i 1900, Schmeltzer [22] published a derivation of the 
covariance of the logarithmic intensity and phase fluctua- . 
tions of a propagating beam of light. His results were 
more general than Tatarski's as they are valid FOr colli - 
mated and focused beams as well as plane and spherical 
waves. However, his results were in the form of compli- 
cated double integral equations and so were not directly 
useable. Fried evaluated these integrals to obtain analytic 
expressions for the logarithmic intensity variance of a spheri- 
cal wave [23] and for a collimated or focused beam [24]. His 
expression for spherical wave scintillation was in agreement 


meen Tatarski's results. 
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mo P969, Ochs [3] repeated Tatarski's experiments in 
Colorado using a helium-neon laser as a light source. His 
experiment found several deviations from Tatarski's theory 
including the existence of a saturation region for the opti- 
cal measurements. When atmospheric turbulence was very 
high, Ochs found that meteorological measurements would 
Sonmeinue to measure turbulence values but the optical 
measurements would saturate at a lower turbulence value. 
Optical and meteorological measurements would not corre- 
fate until pet Srumnewee had decreased below the satura- 
tion level of the optical measurements. 

Similar experiments were performed over land by 
Fitzmaurice in 1970 [4] and by Dowling and Livingston 
in 1971 [2]. Both groups measured the scintillation of 
0.6328 and 10.6 micrometer laser beams and compared optical 
and micrometeorological estimates of atmospheric turbulence. 
Both groups found evidence of a saturation region. 

The work described in this thesis was part of an inter- 
departmental program formed as a combination of the efforts 
of existing optical and meteorological research groups at 
the Naval Postgraduate School. This study applied Tatarski's 
formulation of Rytov's method to describe atmospheric turbu- 
lence. Optical ahd micrometeorological measurements of 
laser beam scintillation and atmospheric temperature fluctu- 
ations have been made independently and then analyzed. Both 


types of measurement yielded values for the index of refraction 
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meructure constant, Cy. Umeery predvets that if the 
turbulence is isotropic and homogeneous, these values of Cy 
should be equal. However, two major unresolved questions 
remained. First, does saturation of optical values of 


C,, exist in the high turbulence regime? If so, what paran- 


N 
eters determine this saturation? Ochs [3] reported that the 
saturation level seemed to vary with mean wind speed and 
the time of day. 

eeeond, Friehe, of Scrzpps Institute, reported recently 
that humidity fluctuations, which were neglected by Tatarski, 
May contribute significantly to refractive index fluctuations 
in the atmosphere [25]. His findings have not been general- 
ly accepted by propagation researchers. The frequency response 
Of available humidiometers is inadequate to check these pre- 
dictions. 

Although considerable evidence exists to support Tatarski's 
propagation theory over land, the problem of propagation in 
the marine environment has not been resolved. If humidity 
Variations influence refractive index fluctuations, laser 
propagation over the ocean may not be described by Tatarski's 
formulation of Rytov's method. This thesis examines the 
validity of Tatarski's propagation theory in the marine 


boundary layer. 
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itt WHeORY 


A. GENERAL REMARKS 

Atmospheric turbulence iS a three-dimensional, nonlinear, 
mrochastic process that occurs in air whenever certain fluid- 
mechanical parameters are exceeded. Random, dissipative 
fluctuations in the homogeneity of the atmosphere will cause 
the velocity profiles of the air to vary in time and space. 
These fluctuations have significant effects on the propaga- 
tion of a laser beam by continuously changing atmospheric 
density. plmts Gauses tne Temractive index to fluctuate 
mtn time and space. 

Atmospheric turbulence (linear) effects on laser beam 
propagation may be classified into two categories. Large- 
scale inhomogeneities will uniformly refract the laser beam 
and cause the beam to wander in space. Small-scale inhomo- 
MemeiciesS Will differentially refract the laser beam and 
cause interference phenomena or scintillation. Small-scale 
turbulence alters the optical path length over different 
segments of the beam and generates interference patterns 
between adjacent wave fronts. This alters the initial phase 
and amplitude distributions of the propagating beam and 
destroys spatial and temporal coherence. 

Due to the inherent causes of turbulence, it is impossible 
Peepredict the transient velocity fields that arise for a 


given set of boundary flow conditions. Statistical analysis 
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must be used to describe atmospheric turbulence. By 
analyzing the variations in atmospheric homogeneity, one 
may describe the fluctuations of the index of refraction 
over a given path. 

Relating the fluctuations in a propagating laser beam 
to atmospheric turbulence is a three-part problem. The 
first part is to find a convenient formulation for charac- 
terizing atmospheric turbulence. This may be done by 
using statistical estimators of the covariance and the 
structure function to describe fluctuations in the index 
@eeretraction. The second part of the problem is to derive 
the equations that govern the propagation of electromagnetic 
waves through a random medium, starting with Maxwell's 
equations. These equations will relate the magnitude and 
Phase of the radiation to the index of refraction of the 
optical path. The final part of the problem is to combine 
the eveistical description of the index of refraction with 
the propagation equation. The result will be an analytic 
relationship Weceinthe Statistical description of the 
refractive index and a Statistical description of the 
phase and amplitude of the propagating beam. 

ime entire derivation is extremely long and difficult. 
Consequently only: the physically most important points of 
the derivation are discussed here, and much of the turbu- 


lence theory and mathematics are omitted. A complete 


treatment of the propagation problem may be found in 


Beeerences 1, 17, 22, 23 and 24. 
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In order to minimize the effects of beam wander on the 
scintillation measurements, the laser beam was diverged. 
The propagating diverged beam was approximated by a spheri- 
cal wave. Thus the general propagation equations obtained 


were evaluated for the case of a spherical wave. 


B. TURBULENCE THEORY 

Wave propagation in a random continuous medium is 
governed by a stochastic wave equation, whose index of 
refraction, n(r,t), is a continuous random point function, 
and characterizes the transmission medium. For line of 
Sight propagation on the earth, optical paths are relatively 
short enabling one to say that the transit time of a given 
Boase x<ront between the source and receiver is much less 
than the time constant associated with changing spatial 
Sestribution of the refractive index (typically ee seconds). 
The time dependence of n(r,t) will be accounted for by 
Metting each Ges ea Je her = (hee. 5... Be conn a 
Single realization of the random process, n(r). Hence, one 
can characterize the transmission medium by a refractive 


index in the form 


Ne 


(1) n(?) = <n(r)) + ny (tr) #1 + ny (4) 


where n, (r) represents the fluctuations in n(r). 
Tatarski [1] assumed that the power spectrum of the 
turbulence, o(K), can be described by the empirical 


Kolomogorov model where 


ZZ 





(2) ® (K) = 0.033 c 2 x ll/s 


2 Z 
N exp(-K /K_ ) 


1s a constant that parameterizes the total amount of 


N 


energy in the turbulence. Tatarski showed that C,, can be 


N 


related to a statistical estimator by 
(3) D(d) = C 


where D (4) is the structure function of the refractive 


index given by 


D_(d) = <{n(r) - n(r+a)]* 


and hence 


(4) ¢,° Cine) fread , 


Assuming isotropic, homogeneous turbulence, the functional 
arguments r and d can be changed from vectors to scalar 
distances. 

The quantity K, the spatial wave number, is given by 
i — 27/2 where & is commonly interpreted, in some sense, 
Mewetne size of the turbulent eddy. This turbulence spectrum 
@S valid for values of K that are in the inertial subrange, 
that is, the range where &y < 2 < Lo. The inner scale, & 


0 
characterizes the size of turbulent eddies where viscous 


0? 


dissipation of turbulent energy becomes predominant (for 
values of & less than La). KH is the spatial wave number 


corresponding to the inner scale (i.e., = 27/2q). Values 


Z5 





meee Creater than the outer scale, Lo» characterize the 
mee2on of the turbulence spectrum where the eddies are not 
homogeneous and where energy is added to the turbulence 
spectrum. For values of £ in the inertial subrange, the 
motion of the larger turbulent eddies causes the formation 
and transfer of energy to smaller eddies with very little 
viscous dissipation. These eddies then form and trans- 
mer energy to smaller eddies and so on until eventually, 


ae approximately & = &., viscous effects become the pre- 


0? 
dominant method of energy dissipation and turbulent energy 

me Gonverted to thermal energy. Typically, for atmospheric 
turbulence in the marine boundary layer, ho Bhs 55 ve epcaleay 
of a few millimeters and Lo 1S OL senle ‘Onder Of a few meters 
(depending on the stability conditions and height above 

the boundary). 


Turbulence theory shows that the covariance of the 


refractive index for homogeneous, stationary conditions 
lm Ci (d) = <n(r+d) n(r)> 


is the three-dimensional Fourier transform of the turbulence 


energy spectrum 
(6) C_(d): = 4n/d fe (K) K sin(Kd) dK . 
O | 
Also, the structure function of the index of refraction, 


given by equations (3) and (4), is related to the turbulence 
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energy spectrum by the transform relation 


2 


Sin (Kd) 
rT” - 


(7) D (4) = Bn ae- )o.(K) K 


under conditions of isotropy, homogeneity and stationarity. 
Given equations (2) and (6), one can relate C,(r) LO Cy. 
This means that Cy can be used to describe the statistics 


of the fluctuations of the refractive index in the atmosphere. 


C. DERIVATION OF THE SCALAR WAVE EQUATION 
The equations describing the propagation of an electro- 


magnetic wave in the atmosphere are derived from Maxwell's 


equations for a non-conducting medium without charge sources, 


(8) (a) Ve D=0 
(b) VeB=0 


(c) VX E= - x= 


- ; . a 
(d) Veg g = oe 


With the additional relations 


(eye D = cE with e = e(r) 


(f) B wtH with 4u 


Uo = constant. 


Taking the curl of equation (8c) one gets 


ant 


: E = nv) xX d 


| 


(9) VXVXE=V(V- ££) - VY 


Q 
ct 


Assuming a time dependence of e°1%t for all field quantities 


and substituting equation (8d) into the right hand side of 
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equation (9), one finds 


2 oy oH) 3 (2D 


(10) V(V : EE) -~ V Ee = 7 ¢ a ne eyes 
_ ae 
— Sie = 
dt 
= eine E 


Starting with Maxwell's equation (8a) and using a well known 


identify for the divergence operator, 


(11) cee Do= Yueh = e¥ - E+ E + Ve = 0 
Or 
(12) en 


Substituting equation (12) into equation (10) gives 


3 ) 2 2 ce = oF e Ve) 
& 


Vv" E + pew 


The index of refraction for a medium is given by 


(14) ae 8 c(ue)+/? = (const) (e) 1/7 


A 
V 


where c equals the velocity of propagation in free space, 
¥V equals the velocity of propagation in the medium and the 


wave number is given by 


(15) 


x" 
M 
N/E 
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mence 
NG 2 

(16) W WE = (=x) (¢ we) = kn and 
Cc 


a7) vee tak eee 


Realizing from equation (14) that i = (const)e, then 
2 


m{ <a 
) 


= = 2Vln(n) and equation (17) becomes 


n 
(18) V° E+ k” n° E = -2V[E += Vin(n)] 


The term -2V[E +» Vin(n)] couples the various components 
of the field together and is responsible for polarization 
Fluctuations in the linearly polarized beam, Fitzmaurice [4] 
Shows that if the deviations of the refractive index from 
its mean value is very small (typically n,/ @) ~ 8S and 
if the wavelength of the propagating wave is much less than 
the inner scale (20) of the turbulence, then one can assume 
that the polarization fluctuation term is negligible. This 


gives the scalar wave equation 
(19) elise Keno E = 0 


where E is any rectangular component: of B(r). 


Za 





me RYTOV'S METHOD 


Pat this point, Rytov [14] transformed equation (19) to 


(20) v*in(E) + \7in(E) |? pen =) 0. > 
Now let 

(21) BE) =r) esplaS(@®)] 

Bnd from equation (1), n(r) = 1 + n, (1). 


Define the complex phase, W, as 

(22) Br) = exp [y(r)} ="exp finA(r) + iS(r}] 
and therefore, 

(23) in(E) = in(A) + iS = w(z). 


Substituting this result into equation (20) gives 


(24) vy + jul? + kG. + n)* = 0. 
Now set 
fe 5) y= Yo a Vy so that 


(26) E 


exp(¥, + ¥,) = E, exp(v,) 


EO is the unscattered wave which satisfies the equation 


(27) a Ea = 0, 
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or the transformed equation 


(28) VIn(E,) + |WIn(E,)|* + k* = 0. 


Substituting equation (22) into equation (28), 
(29) vv, + [Wyyl% + k* = 0 


—_ 


Now, substituting equation (25) into equation (24) gives 


Z 


Z 
| Ps 


a Loe +e? Veet te 


1 = 0 


2 2 : = 
(30) Vp, + Vp, + [YH + VO, 


moe subtracting this from equation (29) gives after some 


manipulation 


+ 2VW_* Vp. + k*(2n + Giese ul) 
1 1 1 


This is the nonlinear Ricatti equation. Since ny is much 
less than one, k’n,’ is much less than 2k*n, and may be 


neglected. Similarly, since [Vv, | << IVvo| ~ k, then 
Jy, 17 is much less. than 2Vy ° Vy, and one can also 


neglect the |p, | term. Finally, one can write equation 


(31) as 


2 = 


(32) vp, + 20y a 


oo VPs sa A'S ny, = 0 


The assumption that Ivy, | << |V),| ~ k is the only 
Significant assumption made in Rytov's method. It is 


equivalent to saying that the change in the variation 


2,9 





of the complex phase, WY,» Over a distance of about one 
wavelength in much less than 2m radians. Realizing that 
the wavelength is much less than the inner scale, dos TG 

is difficult to imagine the phase changing significantly 

in a propagation distance less than hg This approximation 


is generally referred to as the Rytov approximation in the 


literature. 
A new quantity is now defined = 
ue 
(5) W = y,e = YE. 
so that 
| _W _y .-vo 


Substituting equation (34) into equation (32) and expanding 


terms gives 


(35) vow + kW = -2k°n E, 


This differential equation has the solution 
(36) W(t) = f 2k*n,(z') EL (#') GCE - F') av 
vi 


where G(r - r') is the Green's function 


(37) G(r - r') = explik(r - r')) 


- - for @ point s@urce. 
4n(r - r') 
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Equation (38) with n, (r) as a random continuous function 
characterizes the scintillation in the propagating wave if 
n, (r) is known analytically. . 
It is not practical to measure the index of refraction 
at every point over the optical propagation path. Instead 
Mis necessary to statistically characterize n, (r) using 
the covariance function, GaGe). This brings one to the 
third part of the problem. The statistics of v, (r) must 
be calculated starting from equation (38) and the statisti- 
eal description of the covariance of the refractive index 
aS given by equation (6). 
E. APPLICATION OF TURBULENCE THEORY TO THE PROPAGATION 
EQUATION 
After making many transformations, substitutions and 
integrations, Schmeltzer [22] obtained the statistical 
expression for the logarithmic amplitude covariance, c+, (d) 
for a beam propagating in a random medium which Fried [23] 
reduced to 


«2 Z oo 


2 
(39) Contd) = Gr Re { S ds f 6 (KI!) 3, (K 44 =) 


Ks(z - s) 


X [1 - exp ( ee )] dk } 
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fmeaeetne case of a spherical wave. The quantity "z" is 

eveweptical propagation path length,and ''s'' is a dummy 

wareable of integration —— the damensaons of length. 
Fried integrates this expression for a spherical wave 


over s and K to obtain the relation 


f 


Z 


. io ,11/6 C 2 


N where 


0.124 k 


(40)  Cé,(d=0) = C, 


; | 2 2 
(41) Cy, = ([in(A) - 1Inf)]*) 
and A is the amplitude of the wave function, |E(r)]|. 


The logarithmic intensity variance, Cor; 1S a more readily 
measurable quantity 


(42) Co,’ = ([in(1) - nqD}*> 


Where I is the intensity of the radiation (I = WO This , 


(43) Coz’ = ([21n(A) - 21n@ J*) 
= 4 ({n(a) - 1 @)”) 
-4C,,° and 
® Ber = wpoigs pg ®.2?/° 4 
or finally 
(45) Ce ee 


N LI 


The foregoing is a brief summary of the major steps in 


the development of the governing equation for the propagation 


Bd 





Peatistics of a spherical wave. Much of the detailed 
mathematics involved in the derivation has been omitted. 
For the complete derivation, the original papers by 
Megarski [1], Schmedtzer [22], and Fraved [23] should be 
consulted. However, the approximations involved in 
reducing the equations to analytic form are interesting 


and will be analyzed here. 


fe) LHEORETICAL ASSUMPTIONS 

The major assumption that underlies the entire theory 
is that the turbulence energy spectrum can be described by 
the empirical Kolomogorov model [equation (4)]. This model 
was developed in connection with wind tunnel experiments. 
Metecrological measurements by Tatarski [1], by Wyngaard, 
Izumi and Collins [19] and by Panofsky [20] verified that 
it was applicable in the atmospheric boundary layer over 
the inertial subrange.) It was also implicitly assumed that 
only those turbulent eddies in the inertial subrange contribute 
Significantly to beam scintillation. Eddies with 2% > L, would 


0 
tend -to cause beam wander. Eddies with & < L, have little 


0 
effect on the scintillation because viscous forces cause 
rapid collapse. 

The assumption that dX << ho was used to justify ignoring 
the cross polarization term in the scalar wave equation 
(equation 18). The helium-neon laser used in the experiment 
had a wavelength of 0.6328 micrometers which is much less 


than the dimensions of the inner scale (2, was of the order 


of a few millimeters). 
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In order to use the perturbation technique to express 
the variations in the index of refraction, the fluctuations 
in refractive index were assumed to be much smaller than the 
mean value. Meteorological measurements confirm that the 
index of refraction in the atmosphere does not increase 
by more than a few per cent due to density variations; 
hence, this assumption seems justified. 

The assumption that rP2/2o*<<l is a very weak restriction 
at optical wavelengths. For a helium-neon laser and an 
inner scale value of one millimeter, the propagation path 
would have to be greater than approximately 10° Mewers enor 
this assumption to become invalid. 

The final assumption was that turbulence was isotropic 
and homogeneous. This approximation is valid for large 
areas With a smooth regular topography. Ocean conditions 
are conducive to isotropic homogeneous turbulence. This 
approximation means that a point turbulence measurement can 
Characterize the turbulence over an entire propagation path. 

The experimental conditions seemed to satisfy the require- 
eas of Tatarski's, Schmeltzer's and Fried's theory for the 
propagation of a spherical wave in a random medium. A goal 
of this experiment was to test the validity of their propa- 


tion theory in the marine boundary layer. 


G. PROBLEM APPLICATIONS 
The calculation of Cu from equation (45) required values 


of the wavelength of the laser beam, the optical path length 
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moe the Variance of the logarithmic intensity fluctuations 

of the laser beam. The wavelength of the laser was known 
mempbe 0.6328 micrometers. fhe optical propagation path 
Jength was measured on a nautical chart. Log-intensity 
variance was measured using the laser scintillometer which 
mewdescribed in the equipment section of this thesis. The 
calculated values of Cy were designated CxO? LWercating “that 
they were derived from optical measurements. These values of 


C,,. were compared with values of Cy Eiat were calculated trom 


NO 
meteorological measurements (designated Cup) - A brief 
explanation of the theory and equations that describe the 


Saiculation of C 1S appropriate at this point. The 


NT 
meteorological theory is described in detail by Panofsky 


-[20] and by Wyngaard, Izumi and Collins [19]. 


H, METEOROLOGICAL TURBULENCE MEASUREMENTS 

Variations in the index of refraction in the atmosphere 
are caused by variations in the density of the air. These 
variations are due to fluctuations in the temperature, 
pressure, humidity and perhaps gravity TERRE . Empirical 
Meeeurements [4] show that the index of refraction at opti- 


cal wavelengths can be represented by’ the equation 


-6 


(46) n(r) = 1.0 + 77.6 at? + 2: hk 


Me 


) X 10 


where atmospheric pressure, P(r) is in millibars, absolute 


temperature, T(r), is in Kelvins, and wavelength, i, is in 
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micrometers. Evaluation of equation (46) at a wavelength 


of 0.6328 micrometers gives 


P 


(47) (es) S40) = 7.10 Teey X 19°° 


It is generally accepted in atmospheric propagation 
studies that temperature fluctuations are the predominant 


mechanism causing index of refraction fluctuations in the 


atmosphere.* Equation (4) defines ts as 


| 2 
(4) c,2 = Ala) n(r+d) J 


mienwcan also define a temperature structure constant, C." by 


3 2 
2 _ ([{T(r) - T(r+d)]%) 


Rearranging terms gives 





([(79.0 : Te) Gh), a 
(50) , 2 T(r) T(r+D) 


N 2/3 


X 10" °) 37) 
d 


and assuming that the temperature variations are much smaller 


than the mean temperature, then 


= - and 


2 


: Mei htrtd) 9 ir) 


Recently, Friehe of Scripps Institute, in a private 
communication, suggested that humidity variations may contrib- 
ute significantly to refractive index fluctuations under con- 
ditions of high humidity. This has not been generally accepted. 
iMeers difficult to measure humidity fluctuations due to a lack 
of humidiometers that will measure fluctuations with a fre- 
meemey greater than 15 Hz. [25]. 
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((79.0 ee X 107°) [T(r) - T (r+) °\ 
a: s T(r) - TC 


> 2) Cy as SS 
Substituting equation (48) into equation (52) gives 


2 _ P ON eZ 
(53) Cu = 679.0 =I Xora) 


mis the analytic expression for Cur 1s 


IF 6 


(54) Cup = (79.0 =z X ie? ac 


JI, 

In order to calculate Cur using equation (54) measurements 
of mean pressure, mean temperature and Ce as given in equation 
(48) must be made. These measurements will be described in 


the experiment section of this thesis. 
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IV. EXPERIMENTAL APPARATUS AND PROCEDURES 


- 


eres’ section willedescribe the apparatus and procedures 
used to make measurements during the project. The specific 
results and conclusions derived from the measurements will 
be reported in later sections. 

Two major systems were used to measure the data obtained 
mreme theethree main experiments. The scintillometer 
system measured the standard deviation of the logarithmic 


mrtensity distribution, C for the transmitted laser beam. 


sh le 
These measurements furnished the data for the calculation of 
Cro from equation (45). The thermosonde system measured the 


temperature structure constant, C.., mean temperature, T, 


T? 
and the mean pressure, P. These measurements furnished 

the data used in the calculation of Cur from equation (54). 
Meartional experiments were performed to measure the intensity 
meaprlity of the laser, the power spect of the scintilla- 
tion and the degree of truncation of the laser beam by the 
maser Output optics. . These experiments verified the 
scintillometer system's performance and accuracy. The 


equipment used in the experiments will be subsequently 


described in this section. 


A. THE SCINTILLOMETER SYSTEM 
The scintillometer system measured the logarithmic 
intensity variance of the transmitted laser beam. The 


main components of the system are illustrated in Figure l. 
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A low powered helium-neon laser transmitted a beam to an 
Samcal receiver. The receiver detected the incoming beam 
and generated an electrical voltage proportional to the 
instantaneous intensity of the beam. This signal was ana- 
lyzed by the laser scintillometer and recorded on magnetic 
tape. The standard deviation of the logarithmic intensity 
of the laser beam was computed by analogue processing methods. 
This section analyzes in detail the operation of this system. 

The laser output and pointing system controlled the spot 
Size and direction of the beam (Figure 2). The laser beam 
was focused by the first lens to a position near the focal 
point of the output lens. An adjustable mount was pro- 
vided to give the output lens two degrees of translational 
freedom. The spot size was controlled by changing the focal 
point of the output lens with respect to the fixed focal 
point of the first lens. During normal operation, the beam 
was diverged to minimize beam wander effects. A telescope 
with cross hairs, bore-sighted with the laser, provided an 
accurate way of aiming the beam. Beam alignment with the 
emeet was controlled by elevation and azimuthal adjustment 
Serews. A chopper was provided to interrupt the beam for 
ten seconds of every minute during which adjustments were 
made on the background compensator at the receiver. 

After transmission through the ee annen es a portion of 
the scintillating beam plus background radiation was collected 
by the input optics in the detector package (Figure 3). The 


incoming radiation was passed through a 0.63 micrometer 
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bandpass filter and focused by a lens through two small 
apertures. The detector package was aligned in such a 
manner as to pass the laser beam plus background radiation 
through the first aperture to the first silicon photocell. 
The second aperture passed only the background signal which 
was reflected from a small prism to the second photocell. 
The photocells generated an output voltage proportional to 
the intensity of the radiation that impinged upon the photo- 
sensitive surfaces. 

The two signals were then subtracted in a differential 
amplifier (Figure 4) thus cancelling the background signal 
and passing only a signal that was linearly proportional to 
the laser radiation. During operation, the laser beam was 
interrupted for ten seconds every minute, during which time, 
the background signals from the two photocells could be 
Saeetronically adjusted in order to exactly cancel each other. 

The scintillation signal (Figure 5S) was amplified by a 
factor of five and then passed through a logarithmic ampli- 
fier. The logarithmic amplifier computed the logarithm of 
the original signal and sent it to a wide band power ampli- 
fier. The signal level was multiplied by a factor of ten 
and recorded on magnetic tape with an F.M. conversion 
instrumentation tape recorder. The signal was then fed to 
a multi-channel pulse height analyzer which sampled the signal 
level periodically. The pulse amplitude was used as a 
criterion to address one set of 200 analyzer memory channels. 


These channels counted the number of pulses of a given 
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voltage level that were accumulated over a variable time 
period (usually around four minutes). The pulse height 
analyzer output consisted of a CRT display of the number of 
muises of a given Signal level as a function of the logarithm 
of the beam intensity. This corresponded directly to a non- 
normalized, probability distribution of the logarithm of 
the laser beam intensity, ln(1I). The CRT display was 
photographed using a polaroid camera and the photographs 
(example shown in Figure 6) were later analyzed to determine 
Cyo: 

To aid in the analysis, a calibration signal generator 
Seevided direct current calibration signals at 25, 50, 
100 and 200 millivolt levels which were accumulated and 
stored on a second set of 200 pulse height analyzer channels. 
The two sets of 200 channels could be displayed in such 
a manner as to superimpose the calibration signals and the 
Scintillation pulse height distribution. Amplifier and 
pulse height analyzer memory drift were found to be negli- 
gible for the system. 

feesirenents by Tatarski [1], Ochs [3] and Dowling [4] 

have shown that the probability distribution of the logarithm 
of beam intensity is Gaussian. This was expected as a 
direct result of the central limit theorem and the randomness 
of atmospheric turbulence. Realizing that the pulse hemeht 
distribution was Gaussian, then the maximum point corresponded 
to the mean value of the distribution ic the standard devia- 


tion of the logarithmic intensity, C could be measured 


RI? 
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directly from the photographs. Realizing that for a 
Gaussian distribution, the full width at half maximum, X,1/2, 


was related to the standard deviation by 


(56) Coy = 0.425 Xe1/2 


then the standard deviation could be determined from the 


following equation. 


fe5)7.) Cor = 0.425 S(Xp1/2 in cm) 


"Ss" is a scaling parameter that related the distance along 
the abcissa in centimeters to the signal value in log- 
intensity units. S was determined by measuring the distance 
in centimeters between calibration marks and by applying 


the relation 


[1n(V,) - In(V,) J 


X5cm A cm 
V 

~ Th _¢ 
Nh 

EO chi 


where V is the calibration voltage, X is the distance along 
the abcissa in centimeters and AX is the distance in centi- 
meters between two calibration marks, Vy and V,- Once Cy, 
was determined, Cy was computed using equation (45). 
Implicit in the system computation of Cor was the 
approximation that ‘how = (in (1) . This approximation 
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was tested by computing both quantities for typical 
probability distributions of intensity and comparing the 
difference. The error introduced for large values of CxO 
was less than 2% and was correspondingly smaller for lower 
values of Cyo: This error was judged negligible compared 
to the experimental error involved in measuring Cur and 
was ignored. 

The dynamic range and frequency response of the 
scintillometer were wide enough to process all scintilla- 
tion signals without distortion. The logarithmic amplifier 
had a dynamic range of 80 db and distorted voltage levels 
from the differential amplifier that were less than one 
millivolt or greater than ten volts. By choosing an 
appropriate gain setting, the scintillation signal's mean 
value at the differential amplifier output terminal was 
controlled to be approximately 100 mV. All other system 
components had dynamic ranges greater than 80 db. The 
minimum signal levels detected during the experiments were 
around five millivolts while the maximum levels were below 
one volt. System frequency response was flat from dc to 
1.5 KHz. However, the scintillation power spectrum was 
measured to be negligible beyond 300 Hz and hence, the 


Signal frequency response was undistorted. 


B. THE THERMOSONDE SYSTEM 
The meteorological equipment used for measuring Cur 


consisted of a platinum wire thermosonde unit, a barometer, 
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a thermometer and a psychrometer. The platinum wire 
thermosonde unit measured the temperature structure function, 


C., by measuring the changes in the difference of resistance 


‘ag 
between two exposed platinum wires mounted on probes. The 
resistance of each wire was proportional to the temperature 
of the wire. The system electronics used analogue processing 


techniques to compute C., from the equation 


ab 


o 2. {It(r) - TUr+a))4 
, q? 3 


where the value of "d', the probe separation distance, was 
constant. This system had a linear frequency response from 
dc to 1 KHz and a wide dynamic range. The system generated 
a voltage signal corresponding to a which was fed to a 
erip chart recorder. The thermosonde output was recorded 
fea function of time. 

Mean pressure was recorded as a function of time from 
a barometer. Mean temperature was read and recorded from a 
mercury thermometer. Relative humidity was also read and 


recorded from a psychrometer. (C,,. was calculated from 


NT 


measurements of C,,, mean temperature and mean pressure 


Te 
using equation (54). 


C. SCINTILLOMETER PERFORMANCE VERIFICATION EXPERIMENTS 
Three minor experiments were performed in conjunction 
with the scintillation measurements to determine the 


Scintillation power spectrum, laser output intensity 


stability and beam truncation by the output optics. 
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Each experiment was necessary to verify the performance of 
the scintillometer system. 
1. Frequency Response and Dynamic Range 

In order to assure that all laser beam signal levels 
were being accumulated by the pulse height analyzer, the 
dynamic range and frequency response of the system were 
tested. The system response to a wide range of input signal 
levels was found to be logarithmic for signals greater than 
one millivolt and less than ten volts. Typical log-intensity 
distributions ranged from approximately 15 to 400 millivolts 
and all measured distributions were within the system's 
dynamic range. 

Similarly, frequency response tests showed that the 
System would record and process signals with a frequency 
range of dc to 1.5 KHz without distortion. It was then 
necessary to determine the scintillation power spectrum to 
justify the frequency response of the scintillometer system. 
An analogue Fourier spectrum analyzer was designed and 
assembled using a wave analyzer, calibration oscillator and 
an X-Y plotter (Figure 7). 

The FM instrumentation tape recorder fed a measured 
scintillation signal through a power amplifier to a wave 
analyzer. The wave analyzer passed the signal through a 
tunable bandpass filter (bandwidth of 3 Hz) and then 
generated a voltage equal to the rms value of the filtered 


Signal. The rms voltage signal was recorded as a function 
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of time on an X-Y plotter. The system output was a graph 
of the scintillation signal instantaneous relative log- 
intensity at the specific frequency as a function of time 
(Figure 8). The graph was then mechanically integrated 
(using a planimeter) to find the average relative intensity 
Sethe Signal at the specific frequency. By tuning the 
bandpass filter to a variety of frequencies and recording 
the signal intensity, the scintillation inieenonee power 
Spectrum could be measured. 

The tunable bandpass filter had a low frequency 
euc-off at 20 ae The system frequency response was ex- 
memaed to 5 Hz by playing the tape of the scintillation 
Signal at four times the original recording speed and thereby 
increasing the original signal frequency by a factor of four. 

A tunable oscillator was used to calibrate the wave 
analyzer bandpass filter by generating a pure sinusoidal 
wave at the desired frequency and tuning the filter to match 
that frequency. An oscilloscope was used to measure the 
oscillator's output accurately. 

2. Scintillation Spectral Analysis 

The scintillation power spectrum was analyzed to 
verify that the system did not distort the high frequency 
components of the scintillation signal. Spectrum analysis 
of three one-minute duration signals was accomplished using 
the Fourier spectrum analysis system. The three signals 
analyzed represented the maximum, minimum and average values 


of Cy? recorded during the first Monterey Bay experiment. 
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The relative intensities of the various frequency components 
were averaged to give a representative scintillation spectrum 
as shown in Graph VIII. | 

High frequency components above approximately 300 Hz 
maa negligible power levels. The system's high frequency 
distortion limit was approximately 1500 Hz. This was 
determined by the tape recorder frequency response. The 
experiment showed that the high frequency components of the 
scintillation signal were not distorted by the system 
eeeetronics. 

3. Laser Beam Tenet Profile Analysis 

The intensity profile of the laser beam was examined 
to determine the degree of truncation by the output optics. 

A scanning photomultiplier was used to measure beam intensity 
as a function of position in the plane of the output aperture. 
Intensity profiles were measured with and without the output 
lens in position and examined for beam truncation by the 

fens holders. 

The intensity scans showed this profile to be 
essentially Gaussian with negligible truncation. The pro- 
files also indicated that the laser power level varied in 
time, during the warm-up period. Further analysis of 
intensity stability was necessary to show that the variations 
in the radiation received by the scintillometer were due to 
atmospheric propagation effects and not to internal intensity 


Vemmathons. 
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4. Laser Power Stability Analysis 

Laser intensity variations were monitored and 
recorded using a photomultiplier tube and strip chart 
recorder. Measurements of intensity as a function of time 
were taken in the center and wings of the beam profile. 
All measurements showed large fluctuations in beam intensity 
during the warm-up period. These fluctuations were in the 
form of a periodic modulation that gradually decreased in 
amplitude and increased in period. Maximum modulation 
occurred immediately after lasing began. After approxi- 
mately one hour of operation, the intensity stabilized and 
the modulation became negligible. However, it was noticed 
that a 5% modulation could be induced in the beam by blowing 
air across the laser vents after the laser had reached 
thermal equilibrium. The cool air created thermal stresses 
on the laser tube and mirror supports, and caused variations 
Moe the optical gain. | 

Realizing that during the second Monterey Bay experi- 
ment, the laser's operating environment continuously subjected 
ishe ees to thermal stress from meteorological phenomena, 
the validity of the data would be in doubt if the internal 
modulations were large. To test the effects of an outdoor 
environment on thé intensity stability, the laser, photo- 
multiplier and strip chart recorder were Neo en ee enmvhe 
roof of Spanagel Hall and the experiment was repeated. A 
short propagation path minimized atmospheric scintillation 


in the beam, 
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The iZserermtensity “eluctwations stabilized from 35% 
modulation when lasing commenced, to approximately 4% 
modulation after one hour. The level of modulation varied 
with wind velocity but never exceeded 4%. 

Assuming that a 5% modulation was the maximum 
internal intensity variation encountered during the second 
Monterey Bay experiment, the magnitude of the perturbation 


effect on the measured value of C... was calculated from 


NO 
(59) a = (an (1/ (1) 1°) where 
I = J. + AI 
S m 


In the equation, I. is the instantaneous intensity of the 
received beam due to scintillation, and AI was the devia- 
tion from the average intensity due to internal perturba- 
tions of the laser output intensity. Calculations using | 
equation (59) with values of AT = +0.05 (1) and a typical 
preopability distribution of I. Showed that the maximum 


increase in Cu values was 6.4% (computed for aC one 


NO 
: m l/s calculations showed that the perturbation 


cor x 10 
effect would be less than 2.22%. 

These calculations represented the maximum increase 
ial Cro values due to a continuous, random perturbation of 
Al. +0.05 (I) . In actual conditions, the measured average 


intensity perturbation was less than 5%. These estimates 


are an upper bound for possible perturbation effect errors 
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in the data taken on the second and third of May. The 

laser was shielded inside a motel room during the first 
Monterey Bay experiment and perturbation effects should have 
been negligible. Perturbation effects on the laser used 
during the San Nicolas Island experiment by the research 


group from the Naval Missile Center are an unknown quantity. 
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V. DATA ANALYSIS 


The primary goals of the laser propagation research 
peeyect were: 

(1) to develop instrumentation and expertise for 
measurement of laser beam scintillation and micro- 
meteorological turbulence, 

(2) to examine the correlation between micrometeoro- 
logical atmospheric turbulence and laser beam 
scintillation and 

(3) to verify Tatarski's [1] formulation of the theory 
of laser propagation through a turbulent atmosphere, 
in the marine boundary layer. 

Three main experiments were performed between January 
and May of 1973 during which optical and meteorological 
data were collected. The first experiment was conducted 
on January 10. The laser radiation was transmitted across 
Monterey Bay from the Holiday Inn to, a receiver in a labora- 
tory at Hopkins Marine Station (see Figure 9). Micrometeoro- 
logical measurements were made on-board the Naval Postgraduate 
School's Research Vessel Acania. The second experiment was 
conducted at San Nicolas Island on March 3. A laser on the 
Island was aimed at an optical detector mounted on the Acania. 
An Padiitional set of laser scintillation measurements was 
made by a research group from the Naval Missile Center, 
Point Mugu, across Laser Bay on San Nicolas Island (see 
Figure 10). The third experiment was conducted across 
Monterey Bay between the Holiday Inn and Hopkins Marine 


Station on the second and third of.May. Micrometeorological 


measurements were taken on the beach adjacent to the 
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folm@ay Inn. Thesdata collectedduring the three main 
experiments will be reported and analyzed in this section. 
General conclusions based on the results of the three 


experiments will be presented in the following section. 


A. FIRST MONTEREY BAY EXPERIMENT 

ime first Mentercy Bay exper#ment was carried out on 
January 10. The laser output and pointing system was 
placed inside a room in the Monterey Holiday Inn. The beam 
was transmitted a distance of 4.3 kilometers across the Bay 
to the receiver located at Hopkins Marine Station. The 
transmission path was horizontal and approximately Six Mmeegers 
above the surface. This path was closely transverse to the 
prevailing wind direction, and was chosen for this reason 
mma tO Minimize land effects on turbulence homogeneity. The 
meteorological data were taken on-board the R.V. Acania 
which was anchored at approximately the mid point of the 
laser transmission path. The optical and micrometeorological 
Measurements were taken at the same height above the sea 
surface. 


The optical and meteorological estimates of C,, are 


N 
listed in Data Table I and are shown as functions of time on 
Graphs I A and B. Graph II is a scatter plot of the values 

of Cx plotted against Cur: 


The optical measurements of C,, were made in the 


N 
laboratory by processing the tape recorded scintillation 


Signal. The pulse height analyzer accumulated the signal 


a2 











2e =) By se eo 


| ao’; a> > = — => Gere 
7 ae 





= TT 
> a ee 
- = tS 








mmeobability distribution over consecutive five-minute 
Memeods, “This corresponded to,a five-minute integration of 
the laser beam ecdakeaiagdion Signal. Meteorological data 
were recorded on magnetic tape and later transferred to a 
strip chart from which the average values of Cy were 
measured for five-minute periods. Lack of experience in 
operating the scintillometer system caused several inter- 
fuptlons in data collection. 

The experimental error in the optical measurements of 
Cy Was predominantly due to errors involved in measuring the 
moe width of the pulse height distribution. For values of 


8 -1/3 


C,,, between 2 and 5S (X10 “m ), the standard deviation of 


NO 
the measurements was approximately 5%. Instrumentation 
error and the error involved in measuring optical wavelength 
and optical path length were much less than 5%. 

ihe, caperimental «error involwed) um the meteorological 
measurements of Cur was predominantly due to inaccuracies in 
determining the average value of Cy over five-minute intervals 
from strip chart recordings of the instantaneous measurements 
of o- To complicate this process, the thermosonde unit 
had high and low scales for different extremes of turbulence. 
Whenever the scale changed, 30 seconds of data were lost 
due to system transients. The standard deviation was 


estimated as approximately 20% for the C measurements. 


NT 
Errors involved in measuring pressure and temperature and 
internal system errors were negligible compared to averaging 


Saer>. Experimental error bars are indicated at one point 
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on each graph. These: estimates of experimental error 
apply to all points on the graph. 

iimanwiygme the data, 1 1s important to realize the 
fundamental difference between the nature of the meteorologi- 
cal and optical measurements. The meteorological measure- 
ments represent the degree of turbulence at only one point 
along the propagation path. The optical measurements 
represent the integration of turbulence effects on laser 
propagation over the total propagation path. The point 
measurements of turbulence showed greater variations in the 
values of Cur because localized cells of high turbulence 
strongly affected the meteorological sensors. The optical 
measurements did not respond to localized cells of strong 
turbulence unless the dimensions of those cells were sig- 
nificant compared to the propagation path length. The laser 
beam scintillation, in effect, averaged the effects of 
atmospheric turbulence over the entire optical path. This 
1s equivalent to saying that the assumption of homogeneous 
turbulence is not entirely applicable to distances on the 
order One several kilometers. Hence, one should analyze the 


data by comparing the average level of C to the optzcal 


NT 


m 6 
easurements of CxO 


The data showed a strong correlation between the optical 


and meteorological values of C,, with the exception of three 


N 
data regions on Graph I. At approximately 1325, the value 
8 1/3 


of Cy, increased to 12.0 X 10 ~ meters but the corre- 


sponding value of CxO increased to only 4.36 x Toit neers. 
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Mmeain at approximately 1350, the value of Cur increased 


to 15.0 x 107 8m" 1/5 


NO 
esx 10 °m Finally, between 1525 and 1630, the 


whereas C remained steady at 

values of Cur Showed large variations. Corresponding optical 
data remained relatively constant. The average value of 

the meteorological data was approximately equal to the 
optical data. 

may Here been 


These deviations between C and Cy 


NO T 
caused by inhomogeneous turbulence which strongly affected 
the temperature sensors but was not sufficiently wide- 
feread to affect me laser beam propagation.. The fact that 
the average value of Cur (between 1525 and 1630) was approxi- 
mately equal to the relatively steady value of Cx9 tended 
to support this explanation. However, it was also possible 
that the turbulence was homogeneous and that the optical 
values of Cy became saturated. Although the saturation 
effect has been observed in over-land measurements by Ochs 
in 1969 [3], the cause of saturation has not been clearly 
understood or explained. 

Graph II is a scatter plot of the optical and 


meteorological measurements of C The two lines drawn 


NC 
no ont equivalence line delineate the 
equals Cy, within one standard deviation of 


on either side of the C 


region where Cx 


eeorimental error. Sixty-five per cent of the data points 
are Within this region indicating good correlation between 


Measurements of C and Cy The meteorological values 


NO sa 
were syStematically higher than the optical values. This - 
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would be expected in the case of saturation. However, the 
lack of a well-defined saturation level for optical values 
peeciudes a definite conclusion. 

This experiment indicated that the optical and 
meteorological systems that had been developed could 
measure atmospheric turbulence and, for the prevailing 
conditions of the experiment, the local thermal fluctuations 


formed a good predictor of beam behavior. 


B. SAN NICOLAS ISLAND EXPERIMENT 

The next experiment was conducted over the weekend of 
March 3 at San Nicolas Island. Laser beam scintillation 
measurements were taken over a one-kilometer path between 
San Nicolas Island and the R.V. Acania. This experiment 
tested the feasibility of making scintillation measurements 
between the shore and a ship. A helium-neon laser on the 
island was aimed at a shipboard detector package. Meteoro- 
logical measurements were made on-board the Acania. Inde- 
pendent measurements of Cxo were made by the optical 
propagation group from the Naval Missile Center, Point Mugu. 

The data from this experiment are listed in Data Table 
II. Graph III shows the measurement of Cy as a function of 
time. | 

The problem of keeping the laser beam lined up with the 
optical receiver proved very difficult. The laser operator 
aimed the beam at a retro-reflector mounted near the detec- 


tor and used the reflection for sighting corrections. 
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A gimballed mount was provided for the detector which could 
be aligned with the laser beam and thereby reduce the effects 
Sacha pemotien On theedetected beameintensity.. This.system 
proved crude and awkward to operate. By making measurements 
at night, the need for a background compensating photocell 
was eliminated. 

The major problem encountered was that the laser beam 
foulcdenot be pointed at the detector for periods of time 
monger than fave to ten secondsedue to shapemotion. The 
detector was moved around in the beam by the roll and pitch 
of the ship causing a kind of inverse beam wander effect 
and theoretically increasing the measured values of Crrg: 
Variations in laser beam intensity were due to both scintil- 
lation and ship's motion. However, by keeping measurement 
times very short, the effect of ship motion was minimized 
and several log-intensity distributions were accumulated. 
Agreement to within about 40% of optical measurements along 
the beach line, optical measurements from shore to ship and 
meteorological measurements on shipboard lead to some 
Mricidence in the basic measurement techniques. 

This experiment showed, however, that a stabilized 
optical system was necessary to keep the laser beam and 
detector aligned.’ The optical system, as constructed, could 
not reduce ship motion effects to an acceptable lemel. It 
was evident that the system needed to be re-engineered for 


Shipboard use. 
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C. SECOND MONTEREY BAY EXPERIMENT 

The final propagation experiment was conducted on the 
second and third of May across Monterey Bay. As before, 
the laser radiation was transmitted from the Holiday Inn 
to Hopkins Marine Station over a horizontal path, five 
meters above the sea surface. The optical measurements 
were taken as previously described for the first Monterey Bay 
experiment. The meteorological measurements were taken on 
the beach at the Holiday Inn, where the temperature sensors 
were positioned approximately two meters below the optical 
propagation path. The values of Cur were extrapolated for 
the difference in measurement levels. 

1. Data Taken on May 2, 1973 

Administrative and technical problems caused many 

interruptions and delays in data accumulation during the first 
day of the experiment. Numerous beam alignment adjustments 
caused the system to be shut down often for realignment. 
The data collected, especially before 1800, were of ques- 
tionable accuracy. Graphs IV A and B show the optical and 
me eorological measurements of C, plotted as functions of 
time. The average values of the two graphs are approximate- 
ly equal, especially for data taken after 1800. The large 


variations in C are due primarily to beam alignment 


NO 
problems. 
Z. Dawe taken on May 3, 1973 
The data taken on the second day of the experiment 


(Graphs IV C, D, E, F, G, H) showed better correlation 
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between the optical and meteorological estimates of Cy. 


Between 0830 and 1030, the values of Cur increased and 
varied rapidly in time. Corresponding measurements of 
Cx9 Showed a lesser increase and did not reach the high 
values of Cur: This was probably due to the saturation 
ona Cy: Local cells of strong turbulence could have caused 
these effects but it was doubtful that isolated cells would 
have caused sequential increases in a period of two hours. 
Assuming that atmospheric turbulence over the 
propagation path was random, one expects that the number of 
Meealized cells of atmospheric turbulence that affect one 
point on the propagation path during a specific time period 
to be representative of the number of cells that passed 
any point along the path during the same time period. 
imere was no reason to expect that localized turbulence 
would cluster around the thermosonde system. One or two 
Strong increases in Cur dtranpea lenge period of time could 
ee blamed on isolated localized cells of high turbulence 
that would have little effect on laser beam scintillation. 


However, successive, strong increases in C during a 


NT 
fone period of time could only be due to an overall 


increase in atmospheric turbulence. These arguments tend 
po SWpport the saturation explanation for the failure of 


Cx tO Correlate With C 


NT. 


From 1200 to 1600, the values of Cx and Cur 
showed high correlation, while atmospheric turbulence 


remained relatively steady. At approximately 1630, laser. 
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beam alignment was lost and data collection ceased until 
45. During this time period, the values of Cup Showed a 
gradual decrease. When Cuo measurements were re-commenced 
at 1845, values for Cx and Cur Showed a high degree of 
@onrelation. 

3. Data Correlation 


Correlationvof Cx With C was checked by plotting 


NT 


C as a function of C in Graph V (A, B, C) and Graph VI 


NO NT 
(A, B, C, D). Graph V A shows the data measured on the 
second of May. Sixty-one per cent of all the data points 
lie within one standard deviation experimental error of the 
Cu? ene equivalence line. Experimental error bars have been 
plotted on only one point but apply for all points. No 
definite level of saturation was observed. Graphs V, B 
and C, show the data from the second of May sorted according 
to the mean wind speed at the time of measurement. No 
definite trends were noticed. For measurements made when 
the mean wind was less than seven miles per hour, 66% of 
all the data One lie within one standard deviation of 
the Bie Cura equivalence line. For measurements made when 
the mean wind was greater than seven miles per hour, 59% 
of all data points lie within one standard deviation of the 
Cy One cquavalence line. The estimated precision and the 
quality of the data and the lack of.a wide range of wind 
speeds do not justify any significant conclusions concerning 


the effect of mean wind velocity on C,, measurements. 


N 


60 





Erapnevi A 1s a S@aater plot of the Cy measurements 
made on the third of May. A large reduction in technical 
problems led to a higher level Or econil@emec in the optical 
and meteorological data taken on this day. Sixty-three per 
cent of all data points lie within one standard deviation 


er the Cy" equivalence line. Graphs VI B, C, and D 


NT 

are scatter plots of Cy measurements sorted according to the 
Mean Wind velocity at the time of peteoreent For data 
taken during relatively calm periods (Graph VI B), data 
point dispersion is very large. Only 23% of all the points 
lie within one standard deviation experimental error of the 
equivalence line. The data taken during periods when the 
wind speed was greater than three but less than six miles 
per hour (Graph VI, C) show much less dispersion. Seventy- 
nine per cent of the data points are located within one 
Standard deviation experimental error of the Cy? One equiva- 
mence line. This trend toward less dispersion with higher 
wind speed is continued in Graph VI D where data points 
measured during periods when the wind speed was greater 

than six miles per hour are plotted. In this case, 93% 

of all the data points are located within one standard 
deviation experimental error of the equivalence line. These 
data indicate that an increase in mean wind speed was 
accompanied by an increase in the homogeneity of the turbu- 


lence and thus, less dispersion of data points on the 


meatter plots. 
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Vi. °GRNERALD RESUS 


Empirical data from the Monterey Bay experiment indicate 
that the optical and micrometeorological measurements of 
atmospheric turbulence correlated well within experimental 
tolerances. Graph VII shows that 68% of all measurements of 
C and C were within one standard deviation experimental 


NO NT 


error of the C Cur equivalence line. The San Nicolas 


NO 
Island data were excluded from the statistical analysis 
because of the unresolved problem of ship motion. These 
results indicate that it is possible to predict laser beam 
scintillation at a wavelength of 0.63 micrometers in the 
marine boundary layer from measurements of mean pressure, 
mean temperature, and the covariance of the temperature 
Eluvctuations. 

Some evidence indicates that saturation of Cro 
measurements may have occurred for values of Cy greater 
than approximately 5.0 x iiome meters 1/3, However, the 
lack of a definite level of saturation in the scatter plots 
maeeests that definite conclusions on this must await 
additional measurements. 

The correlation between the optical and meteorological 
measurements of the refractive index structure constant 
increased as the mean wind speed increased. Measurements 
taken on the second of May were inconclusive due to a lack 


of range in wind speed values. It was expected that  turbu- 


lence would become more homogeneous as the mean wind increased. 


62 





“a ) ae 





The theoretical model, from which the relations for Cy 
were derived, assumed homogeneous turbulence in the hori- 
Zontel plane. Hence, with an increase in»swind speed, 
turbulence should become more homogeneous and the theoretical 
model should describe the physical processes more accurately 
in accordance with the experimental results. 

The San Nicolas Island experiment demonstrated the 
feasibility of making scintillation measurements between 
the R. V. Acania and a shore station. However, problems 
mlaser alignment due to ship*motion indicated that the 
optical receiver system needed to be redesigned. A 
stabilized control system is needed to keep the laser beam 
and detector system aligned. 

Finally, the high correlation between the optical and 
meteorological measurements of the refractive index structure 
@emstant indicates that Tatarski's description of electro- 
magnetic propagation is essentially valid in the marine 
boundary layer at a wavelength of 0.63 micrometers. The 
malidity of the uHeOTy appears to improve with increasing 
wind speeds. Additional experiments are needed to obtain 
Ponclusive Gvidence to support Tatarski's formulation over 


a wide range of meteorological conditions. 
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VII. ‘FUTURE RESEARCH SUGGESTIONS 


Refractive index structure constant measurements should 
be made over a wide range of turbulence to determine the 
possible existence and description of the saturation of 
optical values. These measurements could further substanti- 
ate or refute the existing results of this project. 

Project measuréments were restricted to a single 
wavelength and propagation distance. Additional experiments 
should be made to determine the validity of the theoretical 
model's predicted dependency on the parameters of wavelength 
and propagation path length. Experiments could be conducted 
at various wavelengths using various types of lasers. 

Equipment should be designed to shicld the laser from 
fluctuating meteorological conditions and to continuously 
monitor the laser's output intensity. This would reduce 
the internal intensity modulation and provide a convenient 
method of determining beam stability. 

Additional thermosonde units should be added to the 
meteorological instrumentation to obtain a more representa- 
tive sampling of turbulence along the propagation path. 

This would reduce the possibility of localized isolated 
cells of inhomogeneous turbulence giving false indications 
of total path turbulence. Additional experiments Should be 
conducted to determine if the shore line boundary has a 


Significant effect on the meteorological measurements. 
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Simultaneous measurements of C could bewmade both on the 


NT 
beach and on-board the Acania and compared for correlation. 

A digital system should be assembled for both optical 
and meteorological data. Digitizing the scintillation 
and temperature fluctuation Signals would greatly reduce the 
experimental error and provide a direct interface with the 
Naval Postgraduate School's IBM computer. 

Finally, the development of a system to make ship-to- 
Shore measurements of laser beam scintillation would be 
highly desirable. The system would require some form of 
meabilized.joptics to reduce ship motion effects. This 


would make the system totally portable and measurements 


could be made almost anywhere on the coast. 
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VIII. “ PROJECT CONCLUSIONS 


The experimental measurements obtained during this 


research project suggest the following conclusions: 


(1) 


(2) 


(3) 


(4) 


(5) 


Systems were developed that made possible simulta- 
neous optical and meteorological measurements of 
the refractive index structure constant in the 
marine boundary layer. 


Laser beam scintillation at 0.63 micrometers was 
predicted from meteorological measurements of 
eemespheric temperature fluctuations, mean 
pressune and méan temperature. 


Tatarski's formulation of electromagnetic 
propagation [1] was valid for propagation in 
the marine boundary layer at a wavelength of 
0.63 micrometers. 


The validity of Tatarski's formulation increases 
with increasing mean wind speed and subsequent 
increasing homogeneity of atmospheric turbulence. 


More research is needed over a wider range of 
meteorological conditions to substantiate the 
generality of the project's conclusions. These 
measurements would resolve the question of the 
possible existence and description of a satura- 
tion effect on optical scintillation measurements 
of the refractive index structure constant. 
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Time 10 January 1973 (PST) 


Time Variation of the Refractive Index Structure Constant 
mS Measured with the Scintillometer System 


Graph I B 
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Time 10 January 1973 (PST) 


Time Variation of the Refractive Index Structure Constant - 
Measured with the Thermosonde System 
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First Monterey Bay Experiment Data Correlation 


Graph II 


Refractive Index Structure Constant 


Measured with the Scintillometer System 
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0.0 6, 6.0 9.0 
Refractive Index Structure Constant Measured with the 
Thermosonde System (X in 


Data taken over Monterey Bay on January 10. Sixty-six per 
cent of all data points are within one standard deviation 
Esperimmental error of the Cx ONT equaumvalence lane. 
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Graph III 





0300 0400 0500 0600 0700 
Time 3 March 1973 (PST) 


Measurements of the refractive index structure constant 
versus time. Data collected during the San Nicolas Island 
Experiment. 


Oo - Cyg measurements by NAVMISCEN optical propagation 


group across Laser Bay on San Nicolas Island. 


(] - Cx measurements made by NPS optical propagation 
group between. the R. V. Acania and San Nicolas Island. 


A - Cyr measurements by NPS meteorological group on- 
board the R. V. Acania. 
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Optical Measurements of Cy versus Time of Day 


Graph IV B 
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Time 2 May 1973 (PST) 


Meteorological Measurements of C,, versus Time of Day. 
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Graph IV C 
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Optical Measurements of Cy versus Time of Day. 
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Optical Measurements of Cy versus Time of Day. 
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Graph VA 
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Optical Measurements of Cy) versus Meteorological Measurements 
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Data taken on the second of May over Monterey Bay. 


Sixty-one "per™cent of"all”data™points are “within one 
standard deviation experimental error of the Cu NT 
equivalence ame, 
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Optical Measurements of Cy versus Meteorological 
Measurements of C.- Data Taken when the Mean Wind 
Speed was less than 7 MPH. 


Data taken on the second of May over Monterey Bay. 
oixty-Ssix per cent of all data points are within 


one standard deviation experimental error of the 
Cy" ONT equiva lew@e. lane. 
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Optical Measurements of C. versus Meteorological 
Measurements of C..Data Taken when the Mean Wind 
Speed was greater than 7 MPH. 


Data taken on the second of May over Monterey Bay. 


Fifty-nine per cent of all data points are within 
one standard deviation experimental emror of thie 
Cu” Cur equivalence line. 
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Graph VI A 
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Optical Medsurements of Cy versus Meteorological 
Measurements of Cy. 


Data taken on the third of May over Monterey Bay. 


omxcty-three per Cent of all data points are within 
one standard deviation experimental error of the 
Cy? Ont equivalence line. 
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Cur (X 10 “m ) 
Optical Measurements of Cy versus Meteorological 
Measurements of C.. Data taken when the mean wind 
speed was less than 3 MPH, 


Data taken over Monterey Bay on the third of May. 


Twenty-three per cent of all data points are within 
one standard deviation experimental error of the 
CQ” Cur equivalence line. 
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Optical Measurements of Cy versus Meteorological 


_ Measurements of Cy. 


Data taken when the mean wind speed was greater 
than 3 MPH but less than 6 MPH. 


Data taken over Monterey Bay on the third of May. 


Seventy-nine per cent of all data points are within 
one standard deviation experimental error of the 
Crug” ONT equivalence line, 
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Optical Measurements of Cy, versus Meteorological 


Measurements of Cy 


Data taken when the mean wind speed was greater than 
6 MPH. 


Data taken over Monterey Bay on the third of May. 


Ninety-three per cent of all data points are within 
one standard deviation experimental error of the 
CyQ7 ONT equivalence fine. 
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Graph VII 
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(x 10787 1/3) 
Optical Measurements eis Cy versus Meteorological 
Measurements of Cy. 
Includes all measurements of Cy made over Monterey 


Bay. Excludes San Nicolas Island measurements. 


wIXty-Cight per cemt of wall data points are within 
one standard deviation experimental error of the 
CxO" Cur equivalence line. 
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LASER OUTPUT AND POINTING SYSTEM 
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APPENDIX C 
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